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We have studied the prospects of measuring the CP violating phase with atmospheric neutrinos 
at a large magnetized iron calorimeter detector considering the muons (directly measurable) of the 
neutrino events generated by a MonteCarlo event generator Nuance. The effect of #13 and 5c p 
appears dominantly neither in atmospheric neutrino oscillation nor in solar neutrino oscillation, but 
appears as subleading in both cases. These are observable in range of E ~ 1 GeV for atmospheric 
neutrino, where solar and atmospheric oscillation couple. In this regime, the quasi-elastic events 
dominate and the energy resolution is very good, but the angular resolution is very poor. Unlike 
beam experiments this poor angular resolution acts against its measurements. However, we find that 
one can be able to distinguish Sep ~ 0° and 180° at 90% confidence level. We find no significant 
sensitivity for Sep ~ 90° or 270°. 
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PACS numbers: 14.60.Pq 

I. INTRODUCTION 

The evidence of neutrino masses and their mixing 

through neutrino oscillation ,lj is the first footprint onto 

the physics beyond the standard model (SM). However, 

the first hint was obtained observing the anomaly by 

1MB in 1986 and then confirmed by Kamiokande in 

1988 0, The oscillation experiments provide the 

measurements of mass squared differences and mixing 

angles. At present 1(3)<7 ranges are 
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0.046). Here 9ij are the mixing angles in Pontecorvo, 
Maki, Nakagawa, Sakata (PMNS) mixing matrix Q and 



to 2 — to 2 . Currently, there is no constraint on 



the CP violating phase Sep or on the sign of Am 31 . The 
origin of CP violation is still an open problem that needs 
both theoretical and experimental exploration. Till now 
the best tests come from only neutral kaon oscillations. 

The mixing angle #13 tells how strongly the atmo- 
spheric and solar oscillations couple and therefore also 
determines the CP violation effects in neutrino oscilla- 
tion. Both #13 and Sep are observable in solar and at- 
mospheric neutrino oscillation experiments only as sub- 
leading. The three flavor effects are masked mainly by 



wide resolutions and systematic uncertainties. There are 
enormous efforts to probe CP violation and to measure 
#13 in long baseline reactor experiments Double Chooz 
@], Daya bay Q, RENO Q and accelerator ex peri ments 
T2K [9] and NO A 0. In a recent paper ll], it is 
shown that there will be only a hint of CP-violation at 
90% confidence level for sin 2 2# 13 > 0.05 for most values 
of Sep if one considers only the upgraded beams for T2K 
and NOvA combined with reactor data. 
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The CP violating phase Sep (Dirac) has been studied 
for a magnetized Iron CALorimeter (ICAL) detector with 
atmospheric neutrinos in but there is no significant 
sensitivity. In this paper, we have re-examined the sensi- 
tivity of Sep- Then we point out the difference from the 
past analysis. The main difference is that they considered 
a high threshold (2 GeV) in their analysis. For events 
beyond this energy, the sensitivity almost vanishes. In 
this work, we consider only muons (directly measurable) 
of the events generated by a MonteCarlo event genera- 
tor Nuance-v3 [3] and the muon threshold energy of 0.8 
GeV. The main advantage of a magnetized detector is 
that it can measure the oscillation effect separately for 
muon neutrinos and anti-neutrinos. This type of detector 
has been proposed by India-based Neutrino Observatory 
(INO) Q. 
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II. ATMOSPHERIC NEUTRINOS 

The atmospheric neutrinos are produced by the in- 
teractions of cosmic rays with the atmosphere. These 
are primarily protons and Heliums and some heavy ions. 
They produce mainly muons, pions and kaons, which de- 
cay into neutrinos. Here we discuss briefly the uncer- 
tainties which arise in flux calculation. For calculation 
of the neutrino flux, one needs the detailed information 
on (i) the primary cosmic-ray spectra at the top of the 
atmosphere, (ii) the hadronic interactions between cos- 
mic rays and atmospheric nuclei, (hi) the propagation 
of cosmic-ray particles inside the atmosphere, and (iv) 
the decay of the secondary particles. The uncertainty 
on primary cosmic-ray spectra has been improved by the 
experiments 
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The hadronic interactions 
are measured from accelerator experiments. However, 
the available data do not cover all required phase space 
for calculation of the neutrino flux. The propagation of 
cosmic-ray particles inside the atmosphere, and the de- 
cay of the secondary particles are handled accurately by 
simulation [l^ . 



Then the total uncertainty can be expressed as 



total 



S 2 K 



where S v is the uncertainty due to uncertainty of tt pro- 
duction in the hadronic interaction model, Sk due to K 
production, 5 a due to hadronic interaction cross section 
and ^ a i r due to atmospheric density profile. 

For calculating atmospheric neutrino flux from pri- 
mary cosmic ray flux, the measured atmospheric muon 
flux is compared with the calculation 19|. The error in 
fi and v fluxes comes mainly in the 7r production in the 
hadronic interaction model. It has been shown that in 
case of pions above 1 GeV fl9| . 

Ac/L A<f> Ve 



However, this equation does not hold for E Js 1 GcV. The 
present estimated uncertainty of atmospheric neutrino 
flux is discussed later in section ITV Al 



III. THE S C p DEPENDENCE OF OSCILLATION 
PROBABILITY 

The algebraic expression for full three flavor oscilla- 
tion probability with CP phase Sep is very long and 
complicated to understand its change with the oscilla- 
tion parameters. We have plotted the oscillogram of 
AP = [P^(S C p = x°) - P^{8 C p = 0°)] in E - L 
plane with x = 90° and 180° for 9 13 = 10°. We set 
|Am§ 2 | = 2.5 x 10~ 3 eV 2 , 6> 23 = 45° and the solar pa- 
rameters at their best- fit values. Since Sep appears to- 
gether with sin in PMNS matrix, the effect Sep in 
P w increases with increase in sin$i3. The solar and at- 
mospheric oscillation are coupled through #13 . The solar 
neutrino oscillation is dominant at E ~ a few tens of 
MeV, while the atmospheric neutrino oscillation is dom- 
inant at E ~ GeV. These are determined from the cor- 
responding mass squared differences. For this reason we 
see that the Sep effect is prominent in sub- GeV region 
and decreases rapidly with increase in energy. Another 
interesting feature of the plots is that over a large regions 
of E — L plane, AP is positive for anti-neutrino and nega- 
tive for neutrino in both cases for inverted hierarchy (IH) 
and normal hierarchy (NH) . This is why the poor resolu- 
tions cannot wash out fully the CP sensitivity. It is also 
expected that the effect of the CP phase is not fully nulli- 
fied due to the systematic uncertainties with pull method 
of chi-square analysis. The reason is following. The sys- 
tematic uncertainties are high compared to the difference 
AP. But, they are constant over a large region of E — L 
plane, while AP changes with L and E. However, from 
these plots one can also find the optimized CP sensitive 
E and L values for the experiments with neutrino beams. 

IV. THE CHI-SQUARE ANALYSIS 

Due to the low statistics at high energy, the \ 2 is cal- 
culated according to the Poisson probability distribution 
defined by the expression: 

n L ,n E r f / n s \ ~1 



X 



E 

2J = 1 



-2N°j In 



\ fe=i / 



fe=i 



The N°, is considered as the number of observed events 
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FIG. 1: The oscillogram of the difference AP = P^Scp = x") ~ P^(S C p = 0°) in E - L plane for x = 90° and 180° with 
#13 = 10° for both neutrino and anti-neutrino with NH. We set other parameter same as fixed for the analysis and given in the text. 



generated by Nuance for a set of oscillation parameters 
with an exposure of 1 Mton.year of ICAL. The Nfj is 
the corresponding number of predicted events (discussed 
later). These are obtained in a 2-dimensional grids in 
the plane of log£ — L QA . The fk is the systematic un- 



certainty of Nfj due to the fcth uncertainty (discussed 
later). The is the pull variable for the fcth systematic 
uncertainty. 

The events are binned in the grid of log E — L QA plane. 
We use total number of log 12 bins tie = 35 (0.8 — 40 
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FIG. 2: The same asQ] but for IH. 



GeV) and the number of L A bins Ul as a function of the 
energy. We consider n L = 2 x 25, 2 x 27, 2 x 29, 2 x 31, 
and 2x33, when E = 0.8-1, 1-2, 2-3, 3-4, and >4 
GeV, respectively. For the down-going events, the bin- 
ning is done by replacing l L 0A ' by ' — L QA \ We choose 
this type of binning to reflect oscillation effect in a bet- 



ter way in the chi-square analysis 20j, |2l| . The binning 



has been optimized in 2CJ, |2l| for the precision study of 



Amg 2 and 023- We consider the mirror L for down going 
events, which is same if the neutrino comes from exactly 
opposite direction. The factor '2' is to consider both up 
and down going cases. We consider number of events 
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in a bin > 4 to maintain x 2 /d-°-f ~ 1 [20f] . We have 
checked that this always keeps it Ss 1.08 at the minima 
in all cases. 

A. Systematic uncertainties 

We divide the systematic uncertainties into two cate- 
gories: I) overall uncertainties (which are independent of 
energy and zenith angle) , and II) tilt uncertainties (which 
are dependent of energy and/or zenith angle). The en- 
ergy dependent flux uncertainty which arises due to the 
uncertainty in spectral indices, can be expressed as 



E 



ME) 



1 + S E log 



E 
10 W r 



b) Similarly, the flux uncertainty as a function of zenith 
angle can be expressed as 



$6 



;„(cos0,) pa $ (cos6y [l + S z \ cos6» zonith - cos 6>g enith |] 

c) overall flux normalization uncertainty 5t N , d) over- 
all neutrino cross section uncertainty 8 a . We consider 
following values of the above systematic uncertainties: 
S E = 5% with E = 1 GcV for E < 1 GeV, and S E = 5% 
with E = 10 GeV for E > 10 GeV, 5 Z = 4% with 

cos6 ,zenith = = 1Q%; and ^ = 15% _ Thege 

derived from the latest calculation of atmospheric neu- 
trino flux 19j . For each set of oscillation parameters we 
minimize the \ 2 with respect to the all pull variables and 
then use these values to calculate the % 2 for that set of 
parameters. 

We consider all uncertainties as a function of recon- 
structed neutrino energy and zenith angle. Here we as- 
sumed that the tilt uncertainties will not be changed too 
much due to the reconstruction. However, on the other 
hand, if any tilt uncertainty arises in reconstructed neu- 
trino events from the reconstruction method, it is then 
accommodated in \ 2 ■ 

To generate the theoretical data for chi-square analy- 
sis, we first generate 500 years un-oscillated data for 1 
Mton detector. From this data we find the energy-angle 
correlated resolutions separately for and (see Fig. 3 
of 2l|) in 35 E v bins (in log scale for the range of 0.8 — 40 
GeV) and 17 cos6^ enith bins (-1 to +1). For a given grid 
of (El — cos^ onith ), we calculate the efficiency of hav- 
ing Eft > 0.8 GeV (threshold of the detector). For each 



set of oscillation parameters, we integrate the oscillated 
atmospheric neutrino flux folding the cross section, ex- 
posure time, target mass, efficiency and resolution func- 
tion to obtain the predicted data in the reconstructed 
\ogE — L 0A grid [26(. We use the Charge Current (CC) 
cross section of Nuance-v3 Q and the Honda flux of 3- 
dimensional scheme This method of obtaining the 

theoretical data in bins of muon energy and zenith an- 
gle has been used and discussed in our previous works 
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We have done this study considering only the muons 
produced in the CC interactions. It is highly expected 
that the addition of hadron energy to the muon energy 
of an event will improve the result. However, since the 
hadron energy resolution depends significantly on the 
thickness of the iron plates, this addition of hadron en- 
ergy will be more reliable in case of GEANT-based stud- 
ies with realistic backgrounds. 

From GEANT simulation of ICAL detector it is found 
that the energy resolution varies 4-10% and angular res- 
olution 4-12% depending on the energy and the direction 
for our considered range of energy. The width of these 
resolutions are very negligible compared to that obtained 
from kinematics of the scattering processes. We have also 
checked from GEANT simulation that the wrong charge 
identification possibility of ICAL detector is also almost 
zero when the magnetic field is ^ 1 Tesla for our consid- 
ered range of muon energy. 

The iron plates are stacked horizontally and the muons 
produced in the horizontal direction or very near to it 
cannot be detected. So we put a selection criteria in our 
analysis. The muons for a given energy must be with 
a zenith angle such that |90° — # zon ith| is greater than 
the half width at half maxima of the scattering angle 
distribution with that energy. This is discussed in detail 
in 0. 



V. RESULTS 

We have marginalized the \ 2 = Xu + xfi over all os- 
cillation parameters Am| 2 , #23, 6*13, <5cp, Am^, and 
6*12 choosing the range of Am| 2 = 2.0 — 3.0 x 10 _3 eV 2 , 
#23 = 37° - 54°, 6»i3 = 0° - 12.5° and 5 C p = 0° - 360°. 
We set the range of Am^ = 7.06 - 8.34 x 10~ 5 eV 2 and 
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FIG. 3: The contours in Sep — 0ia plane for the inputs of IH and 613 = 10° with Sep = 0°, 90° and 180°, respectively. 



012 = 30° - 40°. However, the effect of Am^ comes 
in subleading order in the oscillation probability when 
E ~ GeV. The 2-dimensional 68%, 90%, 99% confidence 
level (CL) allowed parameter spaces (APSs) are obtained 
by considering \ 2 = Xmin + 2-48, 4.83, 9.43, respectively. 
We show the contours in $13 — Sep plane in Fig. [3] for the 
inputs of IH and 6 13 = 10° with 5 C p = 0°, 90° and 180°, 
respectively. We see that both upper and lower bounds 
can be obtained for 5 C p = 180° at 90% CL. In Fig. gj 
we show the contours for #13 = 7.5° with Sep = 180°. 
The contours with inputs of NH and Sep — 180° are also 
shown in Fig. [5J Here, we find that both upper and lower 
bounds can be obtained with atmospheric neutrinos only 
when 5c p ~ 0° and 180°. However, we find no bounds if 
S C p = 90°. 

The absolute bounds on Sep for its different inputs are 
shown in Fig. [SJ This is obtained after marginalization 
over all oscillation parameters except Sep- In Fig. we 
see that the sensitivity falls drastically if one changes the 



threshold of the detector from 0.8 GeV to 2 GeV. This 
result can be understood from Figs. [T]and[2l Now we can 
conclude that the events at energy range J?k1-2 GeV 
contribute mainly in CP phase sensitivity in atmospheric 
neutrino oscillation. 

It should be notable here that the binning of the data 
to obtain the results in Fig. [7] is different from all other 
analyses in this paper. Here, we have used linear binning 
in E and cos6* zon i t h with bin size of 1 GeV and 0.066, 
respectively. However, we have enlarged the bin size of 
E by adding the nearest bin in E if number of events in 
a bin is < 4. This type of binning has been considered 
to compare our results with that obtained in [1 21 ] . If we 
use the same threshold of 2 GeV as used for the analysis 
in . the results are very similar in nature with 12 1. 
However, here we have considered the energy range 2- 
35 GeV, while the range 2-8 GeV is used for analysis in 
[12I ]. From this study, one extra point we gain is that 
results does not change significantly for different types of 



7 




FIG. 5: The same as Fig. [3] but for different inputs as shown in plot. 



binning. However, we have chosen the binning in log-E — 
L QA plane for all other analyses. This was optimized for 
the precision study of Am^ and #23 21]. But, here we 
checked that the type of binning does not change the 
results significantly. 



VI. THE CASE WITH S C p » 90° OR 270° 

At 8c p = 90° or 270°, there is practically no Sep 
dependence in the muon neutrino survival probability 
P(Vfj, — > Vy) in symmetric matter profile [24|. This can be 
understood in the following way. The oscillation proba- 
bilities can be expressed as linear combinations of co- 
sine and sine functions of 6c p 2J,|25|. The probability 
P{ v ii — > in symmetric matter profile can be expressed 



as 



P(v u 



AcosScp + Dcos2Scp + C, 



where, the coefficient A, B and C are functions of oscil- 
lation parameters, but independent of 5c p- The magni- 
tude of D is 0(Arri2i sin 2 #13) and hence it is difficult to 
observe the effect of the term proportional to cos 26c p ■ 



So, effectively, 



is proportional to cos 5c p 



and there is no 5c p dependence when 5c p = 90° or 
270°. The dependency grows very slowly as one goes 
away from these particular values and hence it is very 
difficult to measure 5c p if its true value is around 90° or 
270°. However, due to slow increase of its dependence, 
one can expect the possibility of its measurement at these 
values. In Fig. [31 we have shown the case for 5c p = 90°. 
We see there is almost no sensitivity. 
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FIG. 6: The variation of A\ 2 [x 2 ~ Xmin] with Sep for different inputs of Sep- The \ 2 ' s marginalized over all oscillation parameters 
except Sep- 
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FIG. 7: The variation of A\ 2 [= X ~ Xmin] with Sep keeping other oscillation parameters fixed. We use linear binning in E and 
cos^ zon ith with bin size of 1 GeV and 0.066, respectively for the chi-square analysis of this plot. This binning is similar with that in 
The lower solid curve represents the result with threshold of 2 GeV used in and dashed one with threshold of 0.8 GeV used 
in our analysis. 



We have studied further the possibilities of improving 
the sensitivity in different optimistic ways. First, we try 
to see if any improvement of the systematic uncertainties 
can make it possible to observe the CP violation. For this 
purpose we have changed the uncertainties and see the 
impact of them on APSs. We see from the discussion in 
section QI] that the flux is less known for E < 1 GeV. We 
have assumed a tilt uncertainty Se = 5% with energy for 
this region. We have checked by improving the energy 
tilt uncertainty from 5% to 2% that there is no much 
significant change in CP sensitivity, except a little hump 
around 90°. To check the effect of horizontal/vertical 



flux uncertainty, we have changed the value of S z from 
4% to 2%. Here also, we find no such significant change. 

If we decrease the threshold from 0.8 GeV to 0.6 GeV, 
both CP sensitivity in oscillation probability (see Fig. Q] 
and [2]) and flux increases. But, we have checked that 
there is also no improvement in APSs with decreasing 
the threshold. The fact is that with decrease in energy 
the angular resolution worsens very rapidly and this kills 
the above prospects. 

At Sep ~ 90° or 270°, we see from the discussion at the 
beginning of this section that the CP sensitivity vanishes 
over whole L — E space. It grows gradually as one goes 
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away from these values. In Fig. [T] and [2] this feature is 
reflected and the difference AP appears relatively in the 
lower energy zone compared with 8c p = 180° case. In 
case of atmospheric neutrinos the direction is not fixed. 
The scattering angle between the muon and the neu- 
trino is very large at this low energy and this mainly acts 
against the CP phase measurements for E < 1 GeV. In 
case of the experiments with neutrino beams, the direc- 
tion is known. This gives the main advantage to measure 
the CP phase there. 

VII. CONCLUSION 

In this paper we have studied the sensitivity of a mag- 
netized ICAL detector in measuring the CP phase with 
atmospheric neutrino oscillation. We have presented the 
results for 1 Mton.year exposure of ICAL, which is 10 



years running of the proposed 100 kTon detector. Here, 
we have considered the muons (directly measurable) of 
the events, which are produced by the charge current in- 
teractions generated by neutrino event generator Nuance- 
v3. We performed a marginalized chi-square study con- 
sidering all possible systematic uncertainties. We find 
that one can be able to distinguish Sep ~ 0° and 180° 
at 90% confidence level. However, there is no significant 
sensitivity for 5 C p ~ 90° or S C p ~ 270°. 
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using the energy-angle correlated resolution functions. 



